The Ground and Low Lying Excited States of the Two-Dimensional Hubbard Model 
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We have studied the ground state of the two-dimensional (2D) Hubbard model by using a quantum 
monte method paying special attention to the shell structure effect on finite size clusters. Our 
calculations show there is a gap for spin excitations in the ground state and incommensurate peaks 
at (tt ± 5,n) and (tt,h ±5) in the spin correlation function for a low lying excited state. In the 
ground state, the long range part of the d-wave superconducting correlation function is enhanced 
and the momentum distribution function at the Fermi level (7r,0) is rounded. The gap in spin 
excitations and the momentum distribution function rounding is consistent with opening a d-wave 
superconducting gap in the ground state. 



FAGS numbers; 7I.10.Fd.,74.20.Mn.,75.40.Mg. 

The low temperature properties and ground state of 
the doped two-dimensional (2D) Hubbard model have 
been intensively studied. This activity increased rapidly 
after the discovery of the high- Tc cuprate, because it is 
widely believed that the single band model of the 2D 
Cu02 plane is the essential model for the high- Tc mecha- 
nism. Although some novel scenarios have been proposed 
for the low temperature state and the ground state of the 
doped 2D Hubbard model, little is known with any de- 
gree of certainty. One fact which is exactly known is 
that no superconducting off-diagonal long range order 
(ODLRO) exists at any finite temperature. [Q However 
it is not known if there is superconducting ODLRO in 
the ground state. 

The auxiliary field quantum monte carlo methods have 
been applied to the doped 2D Hubbard model. The re- 
sults obtained so far by using the ground state algorithms 
have not been sufficient to create a consensus whether or 
not there is superconducting ODLRO in the ground state 
of the model. The difficulty partly comes from the 
negative sign problem of quantum monte carlo methods, 
which prohibits one from studying the ground state of 
the doped 2D Hubbard model in wider classes of physi- 
cal parameter ranges, carrier densities, band structures, 
and lattice sizes. Algorithms have been proposed trying 
to reduce the negative sign ratio. Q 

We have studied the ground state of the doped 2D 
Hubbard model by using the recently developed adaptive 
sampling quantum monte carlo (ASQMC) method, with 
which the negative sign ratio can be greatly reduced, j^] 
To simulate the Fermi degeneracy among k points around 
the Fermi level which may be indispensable to enhance 
the d-wave superconducting correlation on the finite size 
lattices, we have carefully tuned the non-interacting band 
structure as well as the carrier density. We found keeping 
the pseudo degeneracy around the Fermi level is quite 
essential to our results. This is a difficult requirement 
for the standard quantum monte carlo methods. 

The Hubbard model we studied has the next- 



nearest hopping term and the third nearest hopping 
term in addition to the nearest neighbor hopping 
term: H = E(jj>(cLcicT + H.C.) + UJ2^n^^n^l - 
t' E(,,,>(4c,.+i^.C.)-t" Ei^,,]Aclc,.+H.C.), where 
the t, t' , and t" terms are the nearest neighbor hopping, 
the next nearest neighbor hopping, and the third nearest 
neighbor hopping terms and the U term is the on-site 
Coulomb repulsion term. We fixed t — 1 and it is the 
unit of energy. Only the Sz = subspace was studied. 

We have studied the Hubbard model on 6 x 6 and 
10 X 10 lattices. Non-interacting band structure and the 
filling are tuned such that the energy gap between the 
lowest unoccupied level (LUL) close to (0.47r, 0.47r) and 
the highest occupied level (HOL) at (tt, 0) : /^lul-hol 
is less than 0.04. The following two cases were stud- 
ied: (i) 6 X 6 lattice with t' = -0.1667, and t" = 0.2. 
The number of electron is 28, (ii) 10 x 10 lattice with 
t' = -0.2, and t" = 0.05. The number of electron is 84. 
In both cases, we introduced anisotropy ±0.0001 on the 
X and y components of t and t" . The band structures 
obtained with these parameter values may correspond to 
those of high- Tc cuprates. |^] It should be noted that 
the fillings do not necessarily correspond to those of real 
high- Tc cuprate and were chosen for the sake of tuning 
^LUL-HOL- We have studied both cases with [/ = 4. 
The U — 6 case was also studied in (i) and some other 
smaller values of U were adopted to check the consistency 
of our results. We used periodic boundary condition in 
the X and y directions. 

In the case of (i) with U — 4, we found there are 
the two quasi-degenerate states with the lowest ener- 
gies. These two states were obtained by using two dif- 
ferent trial wavefunctions; the non-interacting wavefunc- 
tion and spin polarized Unrestricted Hartree-Fock (UHF) 
wavefunction. The energy difference is about 0.01. Holes 
are doped around (tt, 0) and its symmetric points in the 
ground state and they are doped around (O.SStt, O.SStt) 
and its symmetric points in the excited state. The spin 
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correlation function < S{q) >= ^/N J2i j ^^P{~'1 " (^i ~ 
rj)) < Si ■ Sj > has a broad peak at Q = (tt, tt) in 
the ground state but it has incommensurate peaks at 
(tt ± (5, tt) and (tt, tt ± S) in the excited state. 
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FIG. 1. The spin correlation function < S{q) > of the 
ground state (a) and the excited state (b) of 10 x 10 lattice in 
the case of (ii) with [7 = 4. The gap is calculated to be 0.04. 

Similar results were obtained in the case of (ii) with 
U = 4:. The energy gap is calculated to be 0.04. The spin 
correlation functions are shown in Fig. ^ The broad 
peak of spin correlation function at Q suggests strong 
quantum and/or dynamical fluctuation induced by the 
doped hole in the ground state. The incommensurability 
S in the excited state is calculated to be about O.Itt when 
the electron density p = 0.84. The real space spin corre- 
lations < ■ Sj^ > and their magnitude | < • Sj^ > \ in 
the ground and the excited state are plotted in Fig. ||. In 
the ground state, the magnitude of the spin correlation 
decays exponentially, while that of the excited state does 
not. This suggests there is a gap in spin excitations in 
the ground state. 

The momentum distribution function < ti^ > of the 
excited state is qualitatively similar to that of the non- 
interacting model. It has a sharp drop at the Fermi level. 
In the ground state, the drop of < rifc > at the Fermi 
level (tt, 0) and (0, tt) is smeared out and transferred to 
(0.47r, 0.47r) and its symmetric points as shown in Fig. ||. 
< rife > is therefore similar to the momentum distribution 
function in the d-wave superconducting state. ||] 

We have calculated the d-wave superconducting corre- 
lation function: d-wave(i?) =< OgO^ >, where = 

J-(J -J J ] Lfc^' - ) 

%/2^ AT R+si Rl R+sy V2^ AT R+vl Rl R+yv' 
In the ground states of the cases (i) and (ii), we found 

enhancement of d-wave(i?) at large distances due to U. 
d-wave(i2) in the ground state of the case (i) is plotted 
against r in Fig. ^, where R = (3,r). The d-wave cor- 
relation function oi U — 4 is enhanced against that of 
t/ = at the largest distance. We observed the enhance- 
ment is persistent up to J7 = 6. d-wave(i?) in the ground 
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FIG. 2. The magnitude of real space spin correlation func- 
tion I < 5*0 • > I in the ground and the excited states of 
10 X 10 lattice in the case of (ii) with U — 4. The correlation 
function is plotted as a function of r where R = (0, r) in the 
linear-log scale. The inset of this figure is an original plot of 
the spin correlation function < Sg • > as a function of r. 
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FIG. 3. The momentum distribution function < Uk > of 
the ground state (a) and the excited state (b). The system 
and the parameter values are the same as the previous figures. 
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state of the case (ii) with f/ = 4 is plotted against r, 
where R = (4, r) in the same figure. The d-wave corre- 
lation function is enhanced at the two largest distances. 
Therefore the enhancement due to U is persistent against 
increase of lattice size. The gap in spin excitations and 
the reduction of the momentum distribution function at 
(tt, 0) found in the ground state are also consistent with 
a d-wave superconducting ground state. 

We have studied the shell structure and/or the band 
structure effects on d-wave(.R) on the 6x6 and 10 x 10 
lattices. The parameter t' in (i) and (ii) is changed as a 
variable. We put U = 4. The enhancement of d-wave(i?) 
is obtained when t' < —0.16 in the case of (i). The region 
shrinks to -0.2007 < t' < -0.1997 in the case of (ii). 
While the width of the flat region of the non-interacting 
band dispersion near (tt, 0) in the case (i) is about 7r/3, 
that of the case (ii) is almost zero. This may be the 
primary reason why the enhancement is achieved in a 
limited region in the case of (ii). This suggests that the 
stability of the d-wave superconducting state is largely 
influenced by the non-interacting band dispersion near 
(tt, 0) and hence t'/t and t"/t. 
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FIG. 4. The d-wave superconducting correlation function 
d — wave{R) of the ground states of the 6 x 6 (a) and the 
10 X 10 (b) lattices in the case of (i) and (ii) as a function of 
r, where R = {x,r) . U is taken as a variable and denoted 
in the legends of figures, x = 3 in (a) and a; = 4 in (b). 
Open circles, closed squares, and closed triangles are results 
obtained with U = 0, 4, and 6, respectively. 



they become unstable. To restore the phase coherence, 
they form Cooper pairs and become stabilized. This sce- 
nario for superconductivity in the cuprate is consistent 
with our numerical results and has also been discussed 
on the basis of a scaling theory. The reduction of 
< Tik > at (tt, 0) due to U is also consistent with this 
mechanism. 

It should be noted that Alul-hol implicitly used in 
most previous studies, which are often done at closed 
shell fillings, is still too large. Q| Such a situation makes it 
hard for superconductivity to have ODLRO in the ground 
state. Small ^lul-hol were adopted by the present 
author and also by Kuroki and Aoki (KA). ||] KA stud- 
ied the small U/t region {U/t ~ 1). In the small U /t 
region, they observed enhancements of the d-wave super- 
conducting correlation at large distances. 1^ The results 
obtained with larger value of U/t {U/t ~ 4 and 6) here 
may be consistent with their results. 

In conclusion, we have studied the 2D Hubbard model 
with special attention to the shell structure on the fi- 
nite size lattice by using the recently proposed ASQMC 
method. We found two quasi-degenerate states with the 
lowest energies. In the low lying excited state, there are 
incommensurate peaks of the spin correlation function at 
(tt ± i5, tt) and (tt, tt ± 5). In the ground state, there is a 
gap in spin excitations and reduction of the momentum 
distribution function near (tt, 0). We observed enhance- 
ments of the d-wave superconducting correlation function 
in the ground state. The gap in spin excitations and the 
momentum distribution function rounded at the Fermi 
level are consistent with a d-wave ground state. 
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Implications of our results may be summarized as fol- 
lows: (1) The energy difference between the lowest ex- 
cited state with the incommensurate spin correlation and 
the d-wave superconducting ground state is very small. 
The stabilities of these states are changed by the non- 
interacting band structure and the filling. We cannot 
rule out a possibility that these two states may coexist 
in a phase separated region in a high- Tc sample. This 
possibility is also a matter of experimental and theoret- 
ical concern today. j|] (2) In the flat band region near 
(tt, 0), doped holes may lose their phase coherence due 
to U and may not be able to remain itinerate and hence 
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